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Modulated reheating scenario is one of the most attractive models that predict possible detections 
of not only the primordial non-Gaussianity but also the tensor fluctuation through future CMB ob- 
servations such as the Planck satellite, the PolarBeaR and the LiteBIRD satellite experiments. We 
study the baryonic-isocurvature fluctuations in the Affleck-Dine baryogenesis with the modulated re- 
heating scenario. We show that the Affleck-Dine baryogenesis can be consistent with the modulated 
reheating scenario with respect to the current observational constraint on the baryonic-isocurvature 
fluctuations. 

PACS numbers: 98.80.Cq, 98.80.Es 



I. INTRODUCTION 



The observed cosmological perturbation is generated 
from a primordial curvature perturbation. It is known 
that the curvature perturbation can be time-independent 
after the cosmic time approximately becomes one sec- 
ond (t > 1 sec) because the cosmic fluid had been per- 
fectly radiation-dominated until that time [l|. We ex- 
pect the corresponding perturbations can be generated 
from the vacuum fluctuation of a light field, which re- 
quires the mass of the field to be lighter than the Hub- 
ble parameter H during inflation 2]. In original single- 
field slow-roll inflation paradigm, the curvature pertur- 
bation mainly originated from the perturbation of the 
so-called inflaton field which induces the primordial in- 
flation. Then the curvature perturbation had been al- 
ready produced at the initial epoch during inflation and 
is constant. In this deviation from the Gaussian 

statistics of the fluctuation (so called non-Gaussianity) 
is small since the perturbation does not change the rela- 
tion between the energy density p and the pressure den- 
sity P very much through the cosmic history, and the 
non-linearity parameter to express the degree of the non- 
Gaussianity should be the order of the slow-roll param- 
eters /nl ~ e and/or rj ~ 0(1O~ 2 ). Then it means that 
we will not be able to detect the non-Gaussianity. 

When we consider other paradigms by assuming other 
light field <t, the curvature perturbation may originate 
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mainly from the perturbation of a. Although the cur- 
vature perturbation from 8a might be initially negligi- 
ble, it will grow later even after inflation ends because 
it generates a non-adiabatic pressure perturbation. This 
mechanism works at any epochs except for complete ra- 
diation domination or complete matter domination. In 
this case we should be able to detect those different types 
of non-Gaussianity in future. 

So far a variety of models which produce a non- 
adiabatic pressure perturbation and generate non- 
Gaussianity even in terms of slow-roll inflation have 
been reported in this direction; during multi-field infla- 
tion 0, HI , before a second reheating through the curva- 
ton mechanism [fj, at the end of inflation dur- 
ing modulated reheating (liT - flfjj or modulated preheat- 
ing [l(| 53]; a t a modulated phase transition [l8l fl9j ]. 
and a modulated trapping |20[. By observing the non- 
Gaussianity, we will be able to discern a better model 
from other ones. In this paper, especially we focus on 
modulated reheating because of some attractive reasons 
mentioned later. 

The latest constraint on the non-linearity parame- 
ter /nl from the WMAP 7-year data is / NL = 32 ± 
21 (68%C.L.) [2l| for the local type of non-Gaussianity, 
which means a null detection of the non-Gaussianity. On 
the other hand, the Planck satellite is expected to reach 
A/nl ~ 5. |22h24| . Therefore we may be able to detect 
the non-Gaussianity in near future within 5 years. 

Another good indicators to discern a model from oth- 
ers must be the tensor to scalar ratio r for the per- 
turbation. So far WMAP 7-year has reported only its 
upper bound r < 0.2. [2l|, which does not exclude 
even the standard single-field quadratic chaotic inflation 
model. It is expected that the PLANCK satellite will 
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be able to observe r ~ 0.1 [22|, |23[. In addition, it 
should be quite exciting that future ground-based de- 
tectors QUIET+PolarBeaR will reach r ~ 0(0.01) [H, 
and KEK's future CMB satellite experiment, LiteBIRD 
will observe r ~ O(10" 3 ). [H,[2^ 

Considering these current and future situations for the 
detectability of non-Gaussianity and the tensor fluctua- 
tion, for comparison in future it should be necessary to 
keep those known models in proper order with respect 
to the predicted values of r and /nl , respectively. How- 
ever, unfortunately it might be known that there exist 
few models which predict both the large non-Gaussianity 
of /nl ~ O(10) and the large tensor to scalar ratio of 
r ~ O(0. 1). Therefore, it should be required for theo- 
rists in advance to build models which predict both the 
large /nl and the large r enough to be detectable in near 
future. 

For this purpose, the modulated reheating model must 
be attractive because it possibly produces large non- 
Gaussianity and can predict a large tensor to scalar ratio 
if the inflation scale is sufficiently large like chaotic infla- 
tion with its Hubble parameter H ~ 10 13 GeV. 

On the other hand, we also have to check consistencies 
of the modulated reheating scenario with other obser- 
vational constraints. As will be discussed later in de- 
tail, in the modulated reheating scenario the curvature 
perturbation is effectively governed by the fluctuation of 
the reheating temperature after inflation STr/Tr. Since 
most class of viable baryogenesis scenarios in modern cos- 
mology depend on the reheating temperature, the modu- 
lated reheating may induce a large baryonic-isocurvature 
fluctuation. 1 In this paper we consider baryogenesis in 
models with supersymmetric (SUSY) extension of stan- 
dard model, especially so-called Affleck-Dine (AD) mech- 
anism [29L |30| . which is naturally realized even in the 
Minimal Supersymmetric Standard Model (MSSM) 2 and 
agrees with observations in broad parameter regions. 
Since good candidates for the light scalar field a could 
be found in SUSY [ll| or supergravity (the local theory 
of SUSY), this direction of discussion should be natu- 
rally motivated. As we have already raised the question, 
however, it might be nontrivial if the AD baryogenesis 
is consistent with the modulated reheating scenarios be- 
cause the produced baryon number sometimes depends 
on the reheating temperature in the normal parameter 



1 Some class of scenarios for dark-matter production also depend 
on the reheating temperature such as gravitino thermal/non- 
thermal production. In this case the modulated reheating is 
severely constrained by observations of the cold dark matter 
(CDM)-isocurvature fluctuation 1271. 12811 - In this paper we are 
assuming a dark matter such as the lightest neutralino which was 
decoupled from the thermal bath with its appropriate thermal- 
relic density to fit the observation. Then the modulated reheat- 
ing does not produce a sizable CDM-isocurvature fluctuation. 

2 See also Ref. ll'J for another mechanism of baryogenesis in the 
SUSY cosmologies and its compatibility with the modulated sce- 



space. Thus we have a strong motivation to search the al- 
lowed parameter region for the AD baryogenesis to avoid 
the constraint on the baryonic-isocurvature fluctuation 
in the modulated reheating scenario. 

This paper is organized as follows. In Sec. II we show 
the basic picture of the modulated reheating scenarios 
and the conditions for the decay rate of the inflaton field 
where the large non-Gaussianity can be predicted. The 
AD baryogenesis is outlined in Sec. III. In Sec. IV we 
look for the parameter space to agree with the observa- 
tions and also discuss a possible isocurvature fluctuation 
that the AD baryogenesis may originally have. Sec. V is 
devoted to conclusions. 



II. MODULATED REHEATING SCENARIO 

Here, we g ive a brief review of the modulated reheating 
scenario [ll(. In such scenario, we consider the decay 
rate of the inflaton, T, depending on a light scalar field, 
<7, which has a quantum fluctuation during inflation, that 
is, r = r(cr). The e-folding number N = J din a, where 
a is a scale factor, measured between the end of inflation 
at t — ij n f and a time after the end of the complete 
reheating, t c , is given by 



N = In 
= In 



a(t c ) 

o(*inf) 



In 



ajtc) 

a(t re h) 



(1) 



where i re h represents a time at dhia/dt = H = T. 

Let us consider the quadratic inflaton potential, 
V(<f>) oc <fi 2 . In such case, during the inflaton oscillating 
phase after the inflation, the energy density of the Uni- 
verse relying on the inflaton decays as p oc a -3 and the 
Hubble parameter, H, evolves as H oc p 1 / 2 . Since after 
the complete reheating the energy density of the Universe 
is dominated by the radiation (p oc a~ 4 and H oc a~ 2 ), 
the e-folding number given by Eq. (fT]) is rewritten as 
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where we have used H (i rcri ) = T. The fluctuation of a in- 
duces the modulated reheating and hence the fluctuation 
of the e-folding number is given by 
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where r'(cj) = dT(a)/da. In terms of the reheating tem- 
perature Tr oc r 1 / 2 , the above expression can be rewrit- 
ten as 
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Based on 8N formalism [3ll - l37| . the curvature perturba- 
tion on the uniform energy density hypersurface, on 
super-horizon scales is given by the perturbation of the 
e-folding number as £ rs SN and hence we find that in the 
modulated reheating scenario the curvature perturbation 
can be generated due to the fluctuation of the decay rate 
of the inflaton. Up to the second order, we have 



(k5N = -- 
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Hence the power spectrum and the non-linearity param- 
eter /nl defined as 



(C(k)C(k')) = (2ir)^V(k)5W(k + k>) 

C = Clin + g/NLCun i 



are respectively given by 
T(k) -- 



i r\ 2 /H h x 1 



6 r 



linf 
"2^T 



/nl = 5(1-^- 



(G) 
(7) 

(8) 
(9) 



As an example, let us consider an interaction term in 
Lagrangian as C m t D — g{a)(jnl)tl> where g{<j) is a coupling 
constant and tp is a light fermion which interacts with the 
standard particles and its energy finally goes to thermal 
bath. Then we obtain the decay rate of the inflation T 
as 



9M 
8tt 



-7710 



where is the mass of the inflaton. 
dependence of g as 



(10) 



Assuming a a- 



1 + A 



Mr 



(11) 



where A(~ 0(1)) is some constant parameter and M cu t 
is a cutoff scale of the effective theory [l3|, 0, EE EH , the 
non-linearity parameter /nl is given by 



/NL 



-10 2 A 



2 \ - 1 ( a /M cu t 



0.1 



(12) 



where we have assumed that <r/M cut <C 1. The power 
spectrum is also obtained as 



which can be marginally comparable to the observational 
value Vobs(k) — 3x 10 -10 [21| or subdominant for normal 
scales of M cut (~ 10 16 GeV - 10 18 GeV) 3 . 

Hence, in the modulated reheating scenario we can eas- 
ily find that the large non-Gaussianity (/nl ~ O(10)) 
can be generated even for the high energy inflation with 
the large r = 16e~ 0.1(iJ inf /10 13 GeV) 2 when the model 
parameters are appropriately chosen. 



III. AFFLECK-DINE MECHANISM 

AD mechanism [29|, [3(| has been known as one of the 
powerful candidates for the successful baryogenesis mech- 
anism. It can be realized by taking advantage of a flat di- 
rection along which scalar potential vanishes in the global 
SUSY limit. Hereafter we call the complex scalar field 
that parameterizes the flat direction as AD field $ and 
assume that it carries non-zero baryon charge (3. 

Though the scalar potential for the AD field vanishes in 
the global SUSY limit, it is lifted by non-renormalizable 
terms, the SUSY-breaking effect and some other effects. 
Let us consider a non-renormalizable superpotential for 
the AD field given by 



tUnr = 



3)Mf 



(14) 



where M* is the cut-off scale and the positive integer 
n depends on the flat direction. Including the SUSY 
breaking effect, the induced scalar potential reads 



V = V$- 
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2n+4 
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MP 
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(15) 



where is a gravitino mass and as is a complex nu- 
merical factor whose amplitude is of order of unity. V$ 
is the soft SUSY breaking effect that depends on the 
SUSY breaking mechanism. The second term is the F- 
term that comes from non-renormalizable operator W nr . 
The last term represents the interaction between non- 
renormalizable operator and the SUSY breaking sector 
coming from supergravity effect, which breaks the ?7(1) 
baryon symmetry and is called as the A-term. 

During and after inflation, the AD field acquires the 
Hubble induced mass from the interaction between the 
AD field and the inflaton through the supergravity effect, 
which can be negative, 



U H 



(16) 



V{k) 
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2 \2 ( cr /M C ut X 



0.1 



(13) 



3 In the subdominant case, we have to consider a dominant com- 
ponent of the power spectrum of the curvature perturbations, 
for example, the inflaton fluctuation. In such case (so-called 
"mixed" case), the expression of the non-linearity parameter /Wl 
is changed, but we can also expect the large non-Gaussianity |2Sll . 
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where ch is a positive numerical factor of order of unity. 

4 

We must note that even before reheating there can 
be thermal plasma just after the end of inflation as a 
subdominant component of the universe. Its temperature 
can be expressed as [38| 



T : 



(17) 



where Mq is the reduced Planck scale. Since the AD field 
can have interaction with the thermal plasma directly or 
indirectly, it acquires a thermal potential [39], 



thermal 



/i 2 T 2 |$| 2 
a 2 T 4 log 



fe 2 |$| 



(/ii$i>n 



(18) 



At the onset of the oscillation, the AD field acquires an 
angular momentum due to the A-term, which represents 
the baryon asymmetry of the Universe rig , 

n B {t osc ) ~ /3m 3/2 ( J ffoscMr) 2/( " +1) sin(n0 inf + a), (21) 

where 9- m i and a are the phases of $ during inflation and 
the constant as in the third term of Eq. (|15[) . respec- 
tively. Just after the onset of the AD field oscillation, 
s 3 ng is conserved since the CP- violating ^4-term comes 
to ineffective quickly. This is because the AD field value 
continues decreasing with time during the field oscillation 
due to the cosmic expansion. Since the entropy density 
decreases as s a a~ 3 after the reheating if there is no 
late-time entropy production, the baryon-to-entropy ra- 
tio ns/s is conserved. Thus its present value is estimated 



in addition to above terms. Here h is the Yukawa or the 
gauge coupling constant of the AD field, and a g = g 2 /47r 
represents the gauge coupling constant. The upper term 
is the thermal mass from the thermal plasma. On the 
other hand, the lower one (called the thermal logarithmic 
term) represents the two-loop finite temperature effects 
coming from the running of the gauge coupling with the 
non-zero value of the AD field. 5 

Thus, the AD field evolves with the effective potential, 



V c ff = V + Vn + Vthermal- 



(19) 



During and after inflation, when the Hubble parameter 
H is sufficiently large, the AD field settles down to the 
time-dependent potential minimum, 



|$| ~ (HM, 



n\l/(n+l) 



(20) 



and traces its evolution. Note that there can be several 
non-renormalizable operators for the AD field but only 
the one with the smallest n determines the dynamics of 
the AD field. Thus hereafter we consider only smaller n 
(n<3). 

Let us consider the evolution of the AD field further. 
As the Hubble parameter decreases, the Hubble induced 
mass also gets small. Then, when i? 2 sc ~ \V^g\, the AD 
field (more precisely its radial component) starts to oscil- 
late around the origin. Here the dash denotes the deriva- 
tive with respect to <j> = \/2|$|, and hereafter the sub- 
script "osc" indicates that the parameter or the variable 
is evaluated at the onset of the AD field oscillation. 



(-) 

V s Jo 



(3m 3/2 T R 

M 2 G m sc 



(ffoscAC) 2/( " +1) sin(n0 inf + a). 



(22) 

Note that there are four scenarios according to which 
term in the potential drives the AD field oscillation. Thus 
H osc is different from each other [|(|. The first sce- 
nario is that the AD field oscillation is driven by the 
zero-temperature potential V$. In this case, the Hubble 
parameter at the onset of the AD field oscillation is 



H osc ~ m (|$| osc ), 



(23) 



where mo(|$|) = Vo'(|$|). The second one is that it is 
driven by the thermal logarithmic term and 



- Mg+ l T 3(n + l) \ l/(« + 3) 
M 2n I 



(24) 



The third possibility is that it is driven by the thermal 
mass and 



H osc ~ (h MgTr) 



2 a/s 



(25) 



When the Yukawa coupling h is rather small, there is a 
discrepancy in the thermal potential at |$| ~ T/h. Thus, 
if n < 3, it is possible that AD field start to oscillate 
immediately after the AD field value becomes |$| ~ T/h. 
In this case, the Hubble parameter at the onset of the AD 
field oscillation, 



In the case where cjj is negative, the AD field falls into the origin 
quickly and does not affect cosmology. 

The signature of the thermal logarithmic term can be both posi- 
tive or negative depending on the flat direction. For example, the 
LH U flat direction receives a positive thermal logarithmic term 
|4(l . Here we consider only positive thermal logarithmic terms. 
The case with negative thermal logarithmic term is discussed in 
Ref. 14111 in a different context. 



M?, +1 T 



-,2(n+l) \ !/(3-«) 



(26) 



This classification is essential to estimate the baryonic- 
isocurvature perturbation as we will see in the next sec- 
tion. 
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IV. BARYONIC-ISOCURVATURE 
FLUCTUATION 

Let us consider the baryonic-isocurvature fluctuation 
Sb, which is commonly defined as 



Sb = 



5n B 
n B 



6 s 5(jib/s) 
s 



n B /s 



(27) 



In the case where the baryon-to-entropy ratio depends 
on the reheating temperature, the baryonic-isocurvature 
fluctuation can be also generated in the modulated re- 
heating scenario [27| . If we assume ns/s oc T^, we have 



S(n B /s) 5T R 

TiB/S Tr 



(28) 



Hence, in the case where the curvature perturbation orig- 
inates mainly from the modulated reheating mechanism, 
from Eq. (gj and Eq. (JSHJ), we have 



Si 



-3pC 



(29) 



The current observational limit for the anti-correlated 
baryonic-isocurvature fluctuation is roughly given by 
l-WCl < O(0. 1) [H H2] and hence it means that the 
models with p > 0(1) are conflict with current observa- 
tions. From Eq. (|22"|) . we find that in general the present 
baryon-to-entropy ratio generated by the AD mechanism 
strongly depends on the reheating temperature. Thus 
it seems that the AD mechanism and the modulated re- 
heating scenario are incompatible. 

However, when n = 1 and H osc ~ {McT^/M*) 1 / 2 (Eq. 
flZH) or n = 3 and H osc ~ (h 4 M G T^ 3 (h ~ 10~ 5 ) (Eq. 
([2"5])). the dependence of the present baryon-to-entropy 
ratio on the reheating temperature is canceled. If such 
situations are realized, baryonic-isocurvature fluctuation 
is not generated from the modulated reheating. It is non- 
trivial whether there is a parameter space in which the 
Hubble parameter at the onset of the AD field oscillation 
comes to the value above and the baryon-to-entropy ratio 
can be explained the present value, (ns/s)o — 0(1O~ 10 ). 
In fact, there is two sets of parameter spaces in the 
gravity mediated SUSY-breaking mechanism [4fj| : (i) 
n = 1,M* ~ 10 22 - 23 GeV, and T R > 10 7 GeV, and 
(ii) n = 3, M* ~ 10 16 GeV, and T R ~ 10 6 GeV. 6 In the 
case (i) this cut-off scale M* suggests that small light- 
est neutrino mass for the LH U flat direction and in 
the case (ii) it coincides with the GUT scale. In these 
cases, the AD mechanism and the modulated reheating 
scenario are compatible. 

Here we comment on another source of the baryonic- 
isocurvature perturbation. In the absence of the Hubble 
induced A-terui, the baryonic-isocurvature perturbation 



is generated from the fluctuation for the angular compo- 
nent of the AD field [44| . The magnitude of the fluctua- 
tion is given by 



-flinf 

27r|$ in 



_L (Em. 

2tt I M* 



n/(n+l) 



(30) 



where -flinf and |$; n f| are the Hubble parameter and the 
AD field value during inflation. Thus, from Eq. (|22j) and 
Eq. (|27[) . we have the net baryonic-isocurvature pertur- 
bation as 



S B 



2tt 



cot(nt 



inf 
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STr 

T R ' 



(31) 



By considering the observational constraint on the un- 
corrected baryonic isocurvature mode (\Sb/C\ ^ ®0-) 
[2ll l42l | ) originated from the fluctuation for the angular 
component of the AD field (e.g., see Ref. [45|]) we find 
that the Hubble parameter during inflation should be re- 
stricted to be (i) fliaf < 10 13 GeV and (ii) H ild < 10 9 
GeV, respectively. 7 Thus in the case (i) the chaotic in- 
flation with ifinf ~ 10 13 GeV is still allowed by the obser- 
vation in terms of the baryonic-isocurvature fluctuation. 
In such case, we can expect the large tensor-to-scalar ra- 
tio at the detectable level in the future experiments, for 
examples, such as Planck satellite, PolarBeaR or Lite- 
BIRD satellite. 



V. CONCLUSION 

Recently, the primordial non-Gaussianity and the ten- 
sor fluctuations have been focus of attention to provide 
the information about the physics of the early Universe, 
especially, the origin of the observed cosmological pertur- 
bations. There would be a lot of the scenarios generat- 
ing large non-Gaussianity, and the ones generating large 
tensor-to-scalar ratio, respectively at the detectable level 
in the future experiments independently of each other. 

On the other hand, however, the modulated reheating 
scenario is quite attractive since it can produce the large 
non-Gaussianity, and would simultaneously predict the 
large tensor to scalar ratio when the inflation scale were 
large like the case of chaotic inflation. 

Of course, we need to check the consistencies of the 
modulated reheating scenario as a mechanism of gener- 
ating primordial curvature fluctuations with other phe- 
nomena in the early Universe. In this paper, we focused 
on the baryon asymmetry in the Universe. In particular, 
we consider the one of the most promising candidates for 
baryogenesis, the AD mechanism. 

We have shown that the modulated reheating scenario 
is consistent with AD baryogenesis in some sets of model 



In these reheating temperatures, the gravitino problem can be 
avoidable 



7 Here we assume cot(n$i n f + a) ~ 0(1) — 0(0. 1). This is valid 
when there is no fine-tuning for the initial phase of the AD field. 
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parameters, which can be motivated from the physics 
at the high-energy scales, even if the observational con- 
straint for the isocurvature fluctuation and the grav- 
itino problem is imposed on the prediction. This conclu- 
sion is not changed even if we consider the uncorrelated 
baryonic-isocurvaturc mode originated from the fluctua- 
tion of the angular component in the AD field. Therefore 
we will be able to discern this model from others in prin- 
ciple when we detect a large /nl and a large r in near 
future. 
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